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The non-equilibrium properties of a driven quasi-one dimensional superconducting ring subjected to 
a constant electromotive force {emf) is studied. The emf accelerates the superconducting electrons 
until the critical current is reached and a dissipative phase slip occurs that lowers the current. The 
phase slip phenomena is examined as a function of the strength of the emf, thermal noise, and 
normal state resistivity. Numerical and analytic methods are used to make detailed predictions for 
the magnitude of phase slips and subsequent dissipation. 
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I. INTRODUCTION 

When driven away from equilibrium, many systems en- 
counter instabilities leading to new states or phases. Of- 
ten there exists a multiplicity of possible states that can 
be selected near the onset of the instability. The selected 
state may depend on various factors such as the rate at 
which the system is driven through the instability, noise, 
internal excitations, different dissipation mechanisms and 
the system size. 

In this paper the selection of states is studied in driven 
superconducting rings. Many of the phenomena observed 
here are not limited to superconducting rings, but ap- 
pear in many other physical systems ranging from pat- 
tern forming systems @,|jH,| to lasers §. The relative 
simplicity of the superconducting system makes it pos- 
sible to obtain information about some of the general 
questions in driven non-linear systems such as state se- 
lection and the effect of dissipation on the state selection 
process itself. 

The mesoscopic nature of the system, i.e., the super- 
conducting ring has a finite circumference with a finite 
number of accessible states, is fundamental to this prob- 
lem. First, it leads to the existence of a finite number of 
metastable current-carrying states which can compete for 
occupation. It is this competition that lies at the heart of 
the problem. Second, care must be taken to distinguish 
between voltage-driven and current-driven systems. As 
shown by Tarlie et. al. 0, for systems that are not in 
the thermodynamic limit, i.e. mesoscopic systems, the 
choice of ensemble is not free. In this paper we focus 
on voltage-driven systems as opposed to current-driven 
systems. 

In addition to providing a prototype system to study 
various aspects involving driven systems in general, non- 
equilibrium superconductivity is of great interest in its 
own right. Indeed, the current-induced transitions in 



superconducting filaments have been a subject of in- 
tense experimental and theoretical study for almost three 
decades. Ref. || provides a comprehensive review of the 
field. 

We concentrate on the emergence of the dissipative 
phase-slip state [p| Jl^Jll| , p^ , [l3|| in voltage-driven meso- 
scopic systems. When a superconductor (below T c ) is 
driven by a voltage-source, the supercurrent increases un- 
til it reaches a critical value at which point the system 
becomes unstable. Several interesting phenomena may 
then occur: the system will enter the dissipative phase- 
slip state, Joule heating can take place, mode locking, as 
well as other phenomena. Here, the focus is on the onset 
of the instability and its effect on the dynamics of the 
superconducting state. 

The transitions between the current-carrying states 
can take place via two fundamentally different routes: 
i) by a nucleation process involving thermal fluctuations 
and an activation energy barrier, or ii) the system may 
be driven to an instability by an external driving force. 
In the context of nucleation and metastability, the decay 
of persistent currents in thin superconductors is an old 
and extensively studied problem |IJ,|l^,|l6|,[l7j . However, 
the latter |jl3f involves a decay from a point of instability, 
and it is relatively poorly understood. One of the major 
difficulties is this: whereas in the case of nucleation the 
decay is from a metastable state involving thermal acti- 
vation and a saddle point, in the latter case the external 
force drives the system to a point of instability where 
there is no energy barrier left, i.e., the energy landscape 
looks locally flat. In this instance the decay and the fi- 
nal state depend on various factors, such as how fast the 
system was driven, the relative strength of fluctuations, 
internal excitations, and so on. This makes precise theo- 
retical formulation of the problem difficult, since it is not 
possible to use the free energy formulation as in the case 
of metastability [Q . 
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II. THE SYSTEM 

The physical system considered is a quasi-one- 
dimensional superconducting ring of finite circumference, 
i.e., the radius of the cross-section area (S) of the su- 
perconducting filament is much smaller than coherence 
length (£) and magnetic penetration length (A): VS <C 
£(T) and V# < A(T), respectively, see Fig. 0a. When 
the ring is placed in a time-dependent magnetic field, by 
Faraday's law of induction, an emf is induced in the ring. 
From London's equation (E(f) — <9 t [(47rA 2 /c 2 ) J s (f, t)]) 
this leads to a current that increases in time. (Here E 
is the electric field, J s the supercurrent density, and c 
the speed of light.) The time-dependent increase in the 
current cannot continue indefinitely, and eventually the 
current will reach a critical value, at which point the sys- 
tem becomes unstable and a dissipative phase slip will 
occur, resulting in a reduction of the current |l|] by a 
discrete amount. 
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FIG. 1. a) Illustration of a voltage driven superconducting 
ring. The magnetic flux is due to an infinitely long solenoid 
passing through the center of the loop, b) Illustration of the 
current carrying states as uniformly twisted plane waves. At 
a phase slip the amplitude of the helix approaches locally 
zero. The line through the center of the helix represents the 
superconducting wire. 

It is important to reemphasize that the system dy- 
namics in the case under study here, viz. the decay of 
the system from a point of instability, is very different 
from the historically well-studied problem of the decay 
of the system from a point of metastability. The picture 
of the system hopping from one local minimum to the 
next no longer applies. Rather, the picture now is one 
where the system is initially in a locally stable state, but 
as a consequence of the voltage source, the energy land- 
scape evolves in such a way that as the critical current is 
reached, the system finds itself at the top of a hill. When 
this situation is encountered, it is possible that there ex- 
ist a variety of different valleys for the system to flow 
into, each valley leading to a locally stable state. In this 



picture, each of these locally stable states compete for 
occupation. 

To examine this phenomena the Ginzburg-Landau the- 
ory of dirty superconductors will be considered. The 
Ginzburg-Landau free energy functional can be written, 
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where A is the vector potential, 4" is the complex valued 
order parameter, e is the electron charge, m e is the elec- 
tron mass, c is the speed of light, h is Planck's constant, 
and a and b are the expansion coefficients. Since the cur- 
rent is induced in the loop by a time varying magnetic 
flux, the effect of the induced electromotive force (emf) 
must be included in the GL description. By Faraday's 
law of induction, the electrons in the loop are subjected 
to an emf 



(2) 



where £ is the induced emf and <I>(i) is magnetic flux 
through the loop. The magnetic flux and the magnetic 
field are related by 



$ = 



B ■ dSi 



(3) 



where Si is the area of the loop. Eqs. (§) and @ can be 
combined to obtain a relationship between the the vector 
potential and the electric field, i.e., if E x is used to denote 
the tangential component of the field then A x = —E x ct, 
where A x is the tangential component of the vector po- 
tential. This in turn gives A x — —£ct/L, where L is the 
length of the wire. 

The one-dimensional nature of the problem allows sev- 
eral simplifications. First, since the wire is narrow the 
magnetic field generated by the supercurrent does not 
significantly influence the order parameter. This allows 
one to treat the vector potential A x as a parameter in- 
stead of as a dynamical variable. In addition, since the 
magnetic field energy due to the supercurrent is much 
smaller than the energy associated with the order pa- 
rameter, the magnetic field term can be dropped from 
the free energy pl| . Finally, since the radius of the wire 
is less than £ the order parameter is only a function of the 
tangential direction (x). The geometry of the wire im- 
plies periodic boundary conditions, i.e., ^(x) = SSf(x+L). 

For further analysis and computational efficiency it is 
convenient to rewrite the equation in dimensionless form 
using the following transformations: 

= (6/|a|) 1/2 * 
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x' = x/i 

A' = 2e£A/hc 
v' ec EE v ec 2eT GL /h 



(4) 



where £ 2 = Ti 2 / (2m e |a|) and it is implicitly assumed that 
the temperature is below the superconducting transition 
(i.e., a < 0). tql is the Ginzburg-Landau time defined 
as 



TGL 



8k B (T c ~TY 



(5) 



and it is the natural measure for time, i.e., t — > t/rcL- In 
the following, we will work in dimensionless units, i.e., we 
perform the transformations as defined above and drop 
out the primes for convenience. 

The last transformation in Eq. (jij) involves v ec , the 
electrochemical potential generated by the normal cur- 
rent, which will be formally introduced in the next sec- 
tion where the GL theory is extended to include normal 
(Ohmic) current generation. In addition, following the 
scalings in Eqs. ([|), it is natural to measure the length 
in units of the coherence length as £ = L/£. Then the 
rescaled boundary condition reads \&(x) = + £), and 
the dimensionless free energy becomes 



T = f-t* dX { l(dx ~ lAx) " + ^*' 4 } 



(6) 



To describe the dynamics of the superconducting 
condensate, relaxational dynamics are assumed leading 
to the standard stochastic time-dependent Ginzburg- 
Landau (STDGL) equation of motion, i.e., 



~dt 



5JF_ 



(7) 



where rj = r)(x,t) is an uncorrelated Gaussian noise 
source with correlations 

(r)(x,i))=0 
(r)*(x, t)r){x',t')) = 2DS(x - x')8{t ~ t'). 

The angular brackets denote an average, and D is the 
intensity of the noise determined by the fluctuation- 
dissipation j[7) theorem as 



D = 



2irk B T 



(8) 



where H c is the critical field, and H 2 oc (1 — t) 2 , £(T) oc 
(l-t)- 1 / 2 and t = T/T c |j§. 

To make the model numerically more tractable, it 



is convenient to make the transformation Jl7 19 $ — > 

^ e iq(t)x^ where 

q(t) =A X = ul-H, (9) 



where to — ToL^eS/h. This transformation twists, or 
winds, the order parameter along the wire. The effect of 
the transformation is to map the current carrying states 
to twisted plane waves as illustrated in Fig. |]b. After 
the transformation, the periodic boundary condition be- 
comes 

V{£ + x,t) = V{x,t)e iq W, (10) 
and the equation of motion obtained from Eq. (^) reads 



9* _ a 2 * 

dt dx 2 



W - # |*| 2 + U^ujx^ + rj. 



(11) 



This formulation neglects the electrochemical potential 
due to normal current generation at a phase slip center. 
Its inclusion is discussed next. 



A. Electrochemical potential 

Eq. ( |ll| ) would be a sufficient description if the gen- 
eration of a normal current at a phase slip could be 
neglected. This approximation is valid when the nor- 

However, 



mal state resistivity is negligible 13 
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the 



Ginzburg-Landau free energy is only valid for 'dirty' su- 
perconductors in which the normal state resistivity is ap- 
preciable even at low temperatures. One aim of the cur- 
rent study is to examine the effect of the resistive normal 
current to the process. To facilitate this goal the equa- 
tion of motion (i.e., Eq. ([Tl|)) must be generalized to in- 
clude the creation of electrochemical potential gradients 
at phase slip locations. 

A phase slip occurs when the system locally loses su- 
perconductivity and becomes a normal Ohmic conductor. 
As discussed above, below T c the system retains the fully 
superconducting state after making a transition to a state 
of lower current. An important question is the effect of 
the generation of normal current on the dynamics and 
the state selection problem. 

To account for the generation of normal current, the 
time derivative in the STDGL equation of motion must 
be replaced by d/dt + iv ec , where v ec = v ec (x,t) is the 
electrochemical potential generated by the normal cur- 
rent 0^®H>I1- With that substitution, the dimen- 
sionless equation of motion becomes 



d_ 

dt 



* = <9 2 * + * - 1*1 



.tax T , 
i—V + T]. (12) 



Physically, the appearance of the electrochemical po- 
tential is due to local charge imbalance in a superconduc- 
tor. Gorkov p3] was the first to point out that, in a su- 
perconductor the Fermi level, and thus the electrochem- 
ical potential is a local time-dependent variable related 
to the coherence of the superconducting state. Qual- 
itatively, if the local charge balance is disturbed, the 
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Fermi level experiences a local time-dependent pertur- 
bation. This in turn affects the local energy gap. Gorkov 
showed that gauge invariance is preserved, if the order 
parameter depends on time as exp(— 2ijj,Ft/Ti) where hf 
is the Fermi energy. This leads to the second term on the 
left hand side in Eq. @. 

The electrochemical potential can be determined by 
combining charge conservation and Ohm's law in the 
following manner. Charge conservation implies that, 
d x {Jn + Js) = 0, where J„ is the normal current and J s 
is the supercurrent [ p5| . From Ohm's law, i.e., d x v ec = 
—a J n , this can be written 



d 2 v e < 
dx 2 



= a- 



dx ' 



(13) 



where a is a dimensionlcss Ohmic resistivity and can be 
written 



a = Pn/Po 

where p n is the normal state resistivity, 
k B T c (l - t)h 



Po 



■K^{T) 2 e 2 H c {T) 2 ' 



(14) 
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t = T/T c and H C {T) is the critical field. For a dirty 
superconductor |28| this can be written 



p = 0.1455 x 



hp k B T c 



(16) 



U F e 2 H 2 (0) 
where If is the mean free path length. 

III. LINEAR STABILITY ANALYSIS 

The aim of the linear stability analysis is to gain insight 
into the stability of the current carrying state against 
small perturbations, how the perturbations grow or decay 
in time, and how different modes are selected. In general, 
when the total current exceeds the critical supercurrent 
an Eckhaus instability occurs. The Eckhaus instability is 
a longitudinal secondary instability that appears in many 
systems exhibiting spatially periodic patterns ]2(].pT| 

To study the Eckhaus instability in superconducting 
rings the order parameter is linearized around a current 
carrying state by setting ^(x,t) = + S^(x,t), where 
"to = \J\ — q 2 e lqx and q = (u>/£)t. is a current carry- 
ing (or uniformly twisted plane wave) state that is a so- 
lution of Eq. (|l|) in the limit oj/£ < (1 - q 2 c ) 2 /q c « 0.77, 
where q c = 1/V3. This limit is satisfied for the range 
of co/ts considered in this paper (i.e., 2 x 10~ 6 < u>/£ < 
2 x 10 -3 ) . Since the system possesses translational invari- 
ance and admits plane wave solutions, the perturbation 
is given in terms of its Fourier expansion, i.e., 



where a kn (t) is the amplitude of mode n associated with 
wavevector k n — 2irn/£. Substituting into Eq. (|l2]), using 
Eqs. (U) and (|l^) to solve for v ec and linearizing in 5^ 
gives an equation of motion for 5^ or in Fourier space for 
a,k n - Setting Ofc„(i) = a/^e^ 9,0 )* leads to an eigenvalue 
equation which can be solved to give 



Xt(q,a) 



-{l-q 2 ){l + a/2)-kl 
± [(l-q 2 ) 2 (l-a + a 2 /4) 

+Aq 2 (kl + a(l - q 2 )) 1 1 " 



(17) 
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FIG. 2. The eigenvalues as a function q(t) = ut/i for dif- 
ferent a's. The solid line is the n = 1 mode, dotted line n — 2, 
and dashed line n = 3 mode. For small q all the modes are 
stable, i.e., < 0, independent of the value of a. Increasing 
the dissipation (a) increases the growth rate of all phase slips 
modes. The caption in a) shows the boxed area. 

When A„ is negative the corresponding mode is stable, 
fluctuations decay back to zero, and the superconducting 
state persists. When is positive the current carrying 
states are unstable with respect to fluctuations of a finite 
wavevector k n . For the following discussion A~ can be 
neglected as it is negative definite. 

In Fig. ||, A + is shown for the first three modes as a 
function of q for several different values of a. For small q 
all the modes are stable, i.e., A+ < 0. The inset in Fig. ||a 
shows that the modes become unstable sequentially; the 
lowest mode first, then the mode n = 2, and so on. The 
time (t n ) at which a given mode becomes linearly un- 
stable is determined by the condition X n (t n ) = which 
gives 



( 



t n i - -Vl/3 + *£/6. 



(18) 



For a wire of infinite length this time corresponds to the 
time at which the current reaches the critical value, i.e., 
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(w/l)t nl = ^1/3 + fc2/6 ^1/3 and J c = q c (l - 
q 2 c ) = 2/V27. While Eq. (||) implies that single phase 
slip (i.e., n = 1) processes will dominate, this effect is 
offset by the rate of increase of A+, i.e., dX^/dq is an 
increasing function of n. This can be seen in the small a 
limit, i.e., 



, (*£ + 2)(4-A£) 



4 + 5fc2 

or for simplicity in the small k n limit 



a + 0{a 2 )- 



V3(3*£ + 2a + •••). 



(19) 



(20) 



Thus the rate of increase of the positive eigenvalue in- 
creases with n. The situation is somewhat analogous 
to the classic tortoise/hare race if only two modes are 
considered (say n = 1 and 2). In this case the tortoise 
(n = 1) begins the race first since t\ < <2, but the hare 
accelerates faster since d\^/dq\ qi < dX% /dq\ Q2 . To first 
order in a, the effect of dissipation is to increase the 
rate of acceleration of both the tortoise and hare equally. 
Since the tortoise begins the race first this tends to favor 
the tortoise winning the race. In terms of mode analysis, 
increasing the dissipation (i.e., a) increases the probabil- 
ity of a single phase slip (n = 1) occurring over a double 
phase slip (n = 2). 

The linear predictions can be used to estimate the rel- 
ative probabilities of a phase slip of order n occurring. In 
the linear prediction the equal time correlation function 



for the 



where 



mode is 



\a n (t)\ 2 ) = ^-e 2 ^ I dt'e- 2 ^'^ (21) 
* Jo 



a(t,a) = I dt'X+(t',a). 
'0 



(22) 



Following the instability Eq. ( 21 ) describes the evolution 
of the n th mode from the initial current carrying state 
described by ty l n = y/l — q 2 exp[iqx) to the new current 
carrying state described by ^ n = a n exp[i(q— k n )x] where 
a n = yl — (q — k n ) 2 . The quantity 



a n = sj {\a n {t)\ 2 ) / a n 



(23) 



describes the 'distance' from the initial to final n th state 
and can be thought of as an orthogonal coordinate in an 
rt-dimensional space. The unit of measure in this space 
is then 



(24) 



If it is assumed that a phase slip has occurred when d = 1, 
then it is natural to interpret the relative probability of 
an n th order phase slip as 



Pn 



(25) 



Eq. ( p5| ) provides a qualitative picture of the state selec- 
tion process and makes it possible to compare the linear 
theory to numerical results. This will be done in Section 
IV (in particular, see Fig. ||). 

In addition to the dependence of P n on A + , P n also 
depends on the noise strength. While this is not directly 
visible from Eq. ( |25| ) it should be noted that the equation 
d = 1 imposes a D dependence on a n and P n . 

Physically, the noise strength depends on the temper- 
ature of the system via the fluctuation-dissipation theo- 
rem. The intensity of thermal noise increases as T — > T c 
as demonstrated by Eq. (||). Thus, close to T c the rel- 
ative importance of the noise becomes increasingly im- 
portant, whereas away for T c the driving force is domi- 
nant. Since a has no time dependence, the expansion of 
J dtiX^(ti 7 a) leads to the same result as obtained by 
Tar lie and Elder [^3|, i.e., in terms of the intrinsic and 
extrinsic parameters, the instability of order n becomes 
active at time r„ = i(d q X^u>£)~ 1 ^ 2 . 

To summarize, the linear analysis shows that the state 
selection has a subtle dependence on both the applied 
driving force and on the intrinsic properties of the sys- 
tem. It is important to note that this analysis can only 
be expected to give a qualitative description of the pro- 
cess since the analysis does not account for competition 
between the various modes. These results will be com- 
pared with numerical simulations of the stochastic time- 



dependent GL equation in Section IV 



IV. NUMERICAL RESULTS 

The parameters that enter the numerical simulations 
can be estimated by considering typical experimental val- 
ues, such as T c = 3K, T = 0.93T C , H c = 300G, and 
£(0) = y/S = IOOOA. With these values the intensity of 
the noise is D = 10~ 3 , the GL time is tql = 1-4 x 10~ n 
and lu w £/23/iV. In the simulations the tempera- 
ture is fixed and thus the intensity of noise is fixed. 
u) was varied between 0.0001 to 0.1. This corresponds 
to electromotive forces from 2nV to 2/iV. For dirty su- 
perconductors the normal state resistivity can vary be- 
tween 0.01 and 1.0/zf2 — cm and the p varies from 1.0 
to 100.0/iS! — cm. Using these values the dimensionless 
resistivity, a rj 10~ 4 — 1.0, depending on the dimensions 
and the material. 

A simple Euler algorithm was used for the time inte- 
gration of Eq. @, and Eq. (|l|) was solved in Fourier 
space. The complex order parameter was separated into 
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its real and imaginary parts. The simulation parame- 
ters were: L = 64, dx = 0.85, dt = 0.2, where dx and 
dt are the smallest discrete elements of space and time 
respectively. 

A useful parametrization of the length of the system is 
ng = £q c /2TT, where q c = 1/V3 from the Eckhaus analysis 
of the GL equation as discussed above, ng is interpreted 
as the winding number of the order parameter when the 
Eckhaus instability is encountered. For the simulations 
to follow ng — 5 (see Fig. [|). This allows enough com- 
plexity due to interaction between different modes, i.e. 
five modes can compete for occupation, while remaining 
numerically tractable. When computing the probability 
of an n th order phase slip, P n , the averaging was typically 
from 2000 phase slip events (small oj) up to 15,000 phase 
slips (large lo). Simulations performed at large values of 
a and u> (not shown here) often lead to unusual results 
which may be due to numerical inaccuracies. 




b) 



d) 



FIG. 3. A snapshot of a phase slip process with L — 64 and 
ni — 5. a) The current carrying states are uniformly twisted 
plane waves, b) Due to fluctuations, the supercurrent at some 
location along the wire grows slightly faster than in the rest of 
the system, c) At a phase slip the system makes a transition 
back to below the critical current by reducing the number of 
loops in the helix, d) After the phase slip, the system retains 
the perfectly superconducting state with J < J c everywhere, 
z-axis: length, x-y plane: Re^] and Im[^]. 



A. Dynamics of the order parameter 

The dynamics of the order parameter around a phase 
slip is illustrated in Figs. || and [|. Fig. || illustrates that a 
current carrying state is a uniformly twisted plane wave. 
As the current increases the helix becomes more tightly 
wound. Due to fluctuations, there will be weak spots 
where the local supercurrent reaches the critical current 
before the rest of the system. This is the point where the 
amplitude of the order parameter starts to decay rapidly 
toward zero. When |^| 2 — > 0, the phase slip center mo- 
mentarily disconnects the phases to the left and right of 
it, the helix looses a loop, and the supercurrent jumps to 
a lower value. This cycle is repeated periodically. 

In Figs. | and g the behavior described using linear 
analysis in the previous section is clearly visible: as the 
supercurrent increases the absolute value of the order pa- 
rameter, l^l 2 , decreases and at the moment of the phase 
slip approaches zero. After the phase slip the order pa- 
rameter rapidly recovers. Fig. [| demonstrates this be- 
havior. This allows the amplitude to relax toward equi- 
librium in the vicinity of the phase slip center (times ti 
and ^3 in Fig. || and Fig. ||c,d). After a short time the 
wire obtains a uniform current (£4 in Fig. 0). 
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FIG. 4. Order parameter just before and after the phase 
slip. As — > 0, the phase gradient must grow in order to 
maintain a constant current. As the amplitude goes to zero, 
it can 'slip' by a multiple of 2tt and relax to a state of lower 
current. 

To quantify the phase slip events it is useful to con- 
sider two quantities: the spatially averaged supercurrent 
and the winding number, both at time t. The spatially 
averaged supercurrent is given by 



Mt) = 



1 

2f 



J s (^)dx. 



(26) 
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The winding number is a measure of the total phase 
change in the system and can be defined as 



W(t) 



1 

2^ 



dx 



dx 



(27) 



As described above, the order-parameter can change its 
winding number by 2?m, where n = ±1,±2... Only 
changes by an integral multiple of 27r are possible in order 
to preserve the continuity of the order-parameter. This 
also implies that at a single (multiple) phase slip, the 
system removes exactly one (integral multiple) fluxoid. 

Fig. [5] displays the time development of the supercur- 
rent and winding number defined in Eqs. (^6|) and (|2j| ) 
respectively. The electric field drives the current to the 
critical current where an instability occurs and the cur- 
rent jumps to a lower value. As suggested by Fig. [| there 
can be several modes simultaneously present. In the fig- 
ure phase slips of order two dominate but occasionally 
there are jumps of order three. The relative occurrence 
of phase slips of all orders is shown in Fig. || as a function 
of driving force (i.e., u>) for several values of a. 
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FIG. 5. Supercurrent and winding number as a function 
of time. When this figure is compared to the state selection 
probabilities in Fig. (JS[) , it can be seen that the probabilities 
of double and triple phase slips are almost equal, but there is 
still a small probability for single slips. 

As discussed in connection with the linear stability 
analysis the appearance of phase slips of different order 
is a subtle issue. For example, every now and then the 
winding number displays little dips, as if the total phase 
slip was a result of a two stage process. It is instructive to 
look at the state selection probabilities in Fig. |] together 
with the dynamics of the supercurrent and the winding 
number in Fig. |. As seen from Fig. [fj] phase slips of or- 
der n = 1 dominate the process at low driving forces. As 
the driving force is increased phase slips of order n = 2 
become dominant and the shape of the probability curve 
of becomes skewed. Order by order other modes become 



dominant in a similar manner. This is consistent with 
the linear stability analysis as shown in Fig. |[ 
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FIG. 6. State selection probabilities from numerical simu- 
lations as a function of the driving force. The closed circles 
denote single phase slips, open circles doubles, closed squares 
triples, open squares quadruples, and closed triangles phase 
slips of order five. 
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FIG. 7. Upper figure: The winding number as a function of 
time. Lower figure: The square of the amplitude of the order 
parameter taken at four different times inside the boxed area. 
As the winding numbers increases, there are spots where the 
amplitude starts to decay (to). This leads to competition and 
coexistence of several modes, and possible several phase slip 
centers. The parameters correspond to the case when three 
modes, n = 1, 2, 3, are present. 
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The little dips referred to above are a result of compe- 
tition between the modes. As seen in the linear stabil- 
ity analysis, modes of lower order become unstable first 
but the higher order ones grow at a faster rate. This 
leads to competition and crossover effects. This implies 
that the dips in Fig. || are not due to a result of a two- 
stage process where a phase slip of higher order occurs 
via two lower order processes, but instead due to the co- 
existence of different modes with different growth rates. 
Fig. ^ illustrates the complicated nature of the phase slip 
when several modes are simultaneously present. There 
is a competition between different modes, and it is even 
possible for several phase slip centers to exist (almost) 
simultaneously. 
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FIG. 8. Upper figure: The time derivative of |^| 2 at the 
location of the phase slip as a function of time. Lower figure: 
The rate of change the electrochemical potential at the phase 
slip center as a function of time. The simulation parameters 
correspond to a case when single slips dominate almost com- 
pletely. The slice runs from immediately before the phase slip 
to just after it. 

Fig. U shows the time rate of change of |\1/| 2 and the 
electrochemical potential at the phase slip center as a 
function of time, when the mode n = 1 is dominant. 
The time frame is selected in such a way that the figures 
cover the immediate vicinity of the phase slip. At the 
moment of the phase slip, |\&| 2 — 0. After the phase slip, 
\f rapidly recovers its equilibrium value. Since there is 
a constant emf acting on the superconductor, \^\ starts 
to decrease after its recovery. As seen in the lower fig- 
ure, v ec regains its equilibrium value (v ec = 0) at the 
phase slip center considerably slower than <3>. This can 
be seen in the following way. The electrochemical poten- 
tial is zero if the current is uniform throughout the sam- 
ple. However, as seen in Figs. |] and ^, the time required 
to reach a uniform current is much longer than the time 
required for healing of the order parameter at the phase 
slip center. Physically, this corresponds to relaxation of 



the charge imbalance [28 in a superconductor. The re- 
laxation is diffusive &1C,S, with time scales typically 
of order 1(T 9 - l(T 10 s. 

As discussed in the preceding section the electrochem- 
ical potential, or dissipation, changes the probability of 
making an n th order phase slip. This can be seen in 
Fig. H for all values of u). To highlight this feature the 
selection probabilities were numerically estimated as a 
function of a for lu = 1CP 3 and are displayed in Fig. 0. 
In this figure the linear prediction (i.e., Eq. (21-25J)) is 
also included for comparison. While the linear analy- 
sis fails to predict the correct amplitudes for different 
modes, it provides the correct qualitative picture and 
predicts the order in which different modes become dom- 
inant. The quantitative discrepancies stem from two fac- 
tors, first the non- linear terms seem to favor a separation 
between the modes. Fig. ^ shows that once a mode be- 
comes dominant, it quickly suppresses all the other ones. 
However, in linear theory all the allowable modes have 
much higher amplitudes at all values of the driving force 
lu [[50). Second, at the limit a — ► the linear theory pre- 
dicts accurately the crossover points where a new mode 
becomes dominant [^0|, but the presence of dissipation 
(finite a) has a significant effect on that as can be seen 
in Figs. ^ and 0. 
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FIG. 9. State selection probabilities as a function of the 
dissipation (i.e., a) for lu — 1CP 3 . 



B. Power dissipated at a phase slip 

A phase slip is a dissipative process where electrical en- 
ergy is locally converted into heat. This is due to Ohmic 
resistance at the phase slip center. Early experiments 
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p0| , p9]| showed that the differential resistance related to 
the phase slip is temperature independent for a wide 
range of temperatures except very close to T c . This is 
a delicate issue; heating due to Ohmic resistance changes 
the local critical current, and issues related to charge im- 
balance and relaxation may become important M . In the 
following, the heat generated at a phase slip is estimated. 
The normal carriers are assumed to follow Ohm's law. 

The Joule heating law can be used to estimate the heat 
generated at a phase slip. The power generated is 

P = f In ■ EdQ = S 1 j n E x dx, (28) 

J volume J 

where dQ is a volume element, 5* is the cross sectional 
area of the ring, j„ is the dimensional normal current 
density and E x is the electric field along the wire. In 
terms of the electrochemical potential The energy per 
unit volume can then be written 



E 
V 



= £ n 







f 


/ \Jn\ 2 dx 


10 


Jo 



dr' 



(29) 



where, £ = (2H 2 /l)a and J n is the dimensionless normal 
current density. 



per unit mass is 

1 AE 
m AT 

where c is the specific heat. The change in temperature 
is then 







AT = T f 


/ Jldx 


Jo 


Jo 



dr , 



(30) 



where 



2H C 2 

To = a, 

icpm 

p m is the mass density and H 2 oc (1 — t) 2 , where t = T/T c 
p8| . Eq. ( |30"1) can be used to estimate the change in 
temperature due to a phase slip. The linear dependence 
of To on (1 — i) 2 expresses the well known fact |Io| , p9| that 
close to T c the effects of Joule heating are less significant. 
Evaluation of AT requires information about the time 
and the length scales involving v ec , and therefore we have 
not estimated it here. Fig. |Io| shows the accumulated 
energy and power dissipated as a function of time. 
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FIG. 10. Upper figure: Energy (in arbitrary units) dissi- 
pated as a function of time. Lower figure: The corresponding 
power (in arbitrary units) dissipated at phase slips. The fig- 
ures show clearly the quantized nature of dissipation. The 
crossover effects are visible in the width of peaks. The sim- 
ulation parameters were dx = 0.85, dt = 0.2, a = 0.01, 
lu — 0.0064 and D — 10 -3 . These parameters correspond 
to a case when three modes (n = 1,2,3) are active at the 
same time. 

The increase in temperature due to a phase slip can 
be estimated using the heat capacity. The heat capacity 



V. CONCLUSION 

Here, the dynamics of accelerated quasi one- 
dimensional superconductors under the influence of a 
voltage source was studied. A constant emf was used 
to accelerate the supercurrent to the critical current, at 
which point the Eckhaus instability is encountered and 
multiple metastable states can compete for occupation. 
Each of these competing metastable states corresponds 
to a state with a different supercurrent. The transition 
to a new state of lower current involves generation of a 
resistive phase slip center that heals after the phase slip. 
Because the system was driven by a voltage source, it al- 
lowed the study of very general phenomenon, namely the 
relation to the general methods and problems in nonlin- 
ear dynamics, statistical mechanics and pattern forma- 
tion. 

Linear stability analysis was used to investigate the 
Eckhaus instability. It was found out that within the lin- 
ear approximation, the state selection process is a com- 
petition of two factors: the characteristic time at which 
a mode a n (t) becomes unstable, and the growth rates 
of the other modes. For small driving forces, the low- 
order modes have time to grow and dominate the pro- 
cess, whereas for larger driving forces the faster growth 
rates of high-order modes lead to their dominance. In the 
intermediate region the competition leads to crossover. 

Numerical simulations were performed by simulating 
the stochastic time-dependent Ginzburg-Landau equa- 
tion. It was found out that the behavior is consistent 
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with the predictions of the linear analysis. Although the 
behavior was qualitatively similar, nonlinearities and in- 
teraction between the phase slips at higher driving forces 
and higher normal current resistivity lead to differences. 

In spite of the simplicity of the system, it displays rich 
and complex phenomena, and more analytical and nu- 
merical studies are needed. To the authors' knowledge, 
there exists no systematic method to study state selec- 
tion in accelerated systems. Recent work [jl8|,[30|, suggests 
that the path integral method of Onsager and Machlup 
pl| may offer a framework for a systematic study of the 
decay of systems from points of instability when multi- 
ple modes compete for occupation. The extension of this 
work to problems where the dynamical system is evolving 
in time, as is the case here, has not been explored. Addi- 
tionally, future work could explore the two-dimensional 
case numerically. 
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